Exactly-solvable models for atom-molecule hamiltonians 
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We present a family of exactly-solvable generalizations of the Jaynes-Cummings model involving 
the interaction of an ensemble of SU(2) or SU(1,1) quasi-spins with a single boson field. They are 
obtained from the trigonometric Richardson-Gaudin models by replacing one of the SU(2) or SU(1,1) 
degrees of freedom by an ideal boson. Application to a system of bosonic atoms and molecules is 
reported. 
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The Jaynes-Cummings (JC) model ,lj provides a sim- 
ple description of the interaction of matter with a radi- 
ation field. It treats a two-level atom in terms of the 
spin-1/2 generators of the SU(2) algebra and describes 
its coupling to a single radiation field in the so-called 
Rotating Wave Approximation. Despite its simplicity, 
the model has had enormous success in quantum optics, 
finding realization in experiments with Rydbcrg atoms 
in microwave cavities @ and optical cavities Q • 

There have also been several extensions of the JC 
model that have likewise proven useful. One example 
is the Tavis-Cummings (TC) model 4], in which the 
spin-1/2 operators are replaced by operators for arbi- 
trary spin, permitting the description of a collection of 
equivalent two-level atoms with a radiation field. This 
model has been solved exactly using the Quantum In- 
verse Scattering Method for and arbitrary SU(2) spin pj. 
Another example is the Buck-Sukumar (BS) model [3, in 
which a specific non-linear interaction between the atoms 
and the radiation field is included and which is also ex- 
actly solvable. More general non-linear terms have also 
been discussed, but they can only be treated approxi- 
mately 8] . A third example is an exactly-solvable atom- 
molecule hamiltonian that describes the photoassociation 
of pairs of condensed bosonic atoms - based on the alge- 
bra SU(1,1), rather than SU(2) - into molecules with a 
linear interaction [3 , or with a non-linear interaction [6| . 

Similar physics is also at play when a molecular Bose- 
Einstcin condensate is produced through photoassocia- 
tion (with or without the interaction with a Feshbach 
resonance) in dilute fermion [Ij| or boson [ll| gases. The 
production of degenerate bosonic Sodium atom-molecule 
mixtures has been recently reported [l2( , though it is still 
an open question as to whether the molecules formed a 
BEC. Mixtures of fermion atoms and molecular dimers 
are better candidates for constructing a molecular BEC 
due to the suppression of molecular decay by Pauli block- 
ing. Indeed, two groups have reported the manufacture 
of molecular BECs from 40 Ka [Rj and 6 Li [lj| fermionic 
atoms, respectively. Unfortunately, no current exactly- 
solvable model can describe these physical processes. 

In this letter, we show how to generalize the JC model 



to accommodate these physical scenarios, as well as oth- 
ers, in the context of exactly-solvable models. The exten- 
sion is to a family of models that involve an ensemble of 
SU(2) or SU(1,1) quasi-spins and a single bosonic mode. 
For the purposes of our discussion, the SU(2) models de- 
scribe fermion pairs and their coupling to a bosonic mode, 
whereas the SU(1,1) models describe the corresponding 
physics of bosonic pairs. The SU(2) models could also 
be used to describe two-level atoms and a bosonic mode, 
but we will not discuss such models here. The general- 
izations we will describe build on the recently-proposed 
Richardson-Gaudin (RG) integrable models [l5j (for a re- 
cent review see 16]). Following the presentation of the 
models, we will discuss their specific application to a mix- 
ture of bosonic atoms and molecular dimers. 

We begin by introducing the generators of the SU(2) 
and SU(1,1) algebras, Kf, K+ and Kr = {K+)\ which 
satisfy the commutation relations 
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[K+,Kr]= T 25 lj K° 



(1) 



The upper sign refers to the bosonic SU(1,1) algebra and 
the lower sign to the fermionic SU(2) algebra, as they 
will throughout this presentation. 

In the quasi-spin or pair representation of the SU(2) 
and SU(1,1) algebras, the generators are realized in terms 
of particle creation and annihilation operators as 
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m ( a i™) creates (annihilates) a boson or a fermion 
in the state \jm), \jfn) is the state obtained by acting 
with the time reversal operator on \jm), and flj is the 
total degeneracy of single-particle level j. 

There are three families of fully integrable and exactly- 
solvable RG models that derive from the SU(2) and 
SU(1,1) algebras, the rational, trigonometric and hyper- 
bolic models, respectively [l5|. For all three families, it 



has been shown how to write the complete set of com- 
muting operators (the integrals of motion) and the cor- 
responding eigenvalues and eigenvectors. Here we focus 
on the trigonometric family, for which the integrals of 
motion can be expressed in terms of the generators as 



Ri = Kf + 2gJ2i 
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For each degree of freedom i, there is one real arbitrary 
parameter r]i that enters the integrals of motion. These 
operators commute amongst themselves and each com- 
mutes with the conserved quantity K° = J^ Kf, which 
is related to the total number operator. 

We now consider the eigenvalue equation for the in- 
tegrals of motion, Ri\$f) = r^), in the seniority-zero 
sector, namely when all particles are paired. Solutions 
with broken pairs can also be readily obtained, as in |15| . 

In this sector, the eigenstates of Ri are given by 
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where |0) is a state that is annihilated by all the K^ and 
M is the number of pairs. The structure of the collective 
operators B\ is determined by a set of M parameters e Q , 
which satisfy the set of coupled nonlinear equations 
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The associated eigenvalues take the form 
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The important point to note here is that any hamilto- 
nian that can be written solely in terms of the integrals 
of motion Ri is likewise exactly solvable, with precisely 
the same eigenvectors as in (0J and with eigenvalues that 
are obtained directly from those in JBJ). 

We now discuss how to construct an appropriate sub- 
set of the trigonometric RG models that are of relevance 
to the quantum problems described in the introduction, 
involving the interplay of a set of SU(2) or SU(1,1) sys- 
tems with a single bosonic mode. To do this, we use 
a trick originally proposed by Gaudin [17J], which in- 
volves replacing one SU(2) or SU(1,1) degree of freedom 
by an ideal boson. For specificity, we denote the SU(2) 
or SU(1,1) degree of freedom to be bosonized as i = 



and the remaining as i = 1, . . . , L. In the limit Slo - 
the generators map onto ideal bosons according to 
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We now introduce a change of notation for the trigono- 
metric functions that appear in |J3J| for the selected degree 
of freedom, Wj = 1/ sin(j]Q — rjj) , Vj — cot(ryo — Vj) > with 
wf — vf = 1. Moreover, we expand these amplitudes in 
the inverse of the divergent degeneracy Qq , 
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thereby introducing a new set of parameters e; (I = 
1,...,L) to replace the r/i 's. 

Inserting J7J and (JHJ) into ©, we obtain new integrals 
of motion that involve the ideal boson degree of freedom, 
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where G = g 

After bosonization of the selected degree of freedom, 
the resulting Ri still satisfy the conditions for an inte- 
grable model. They remain hermitian, global, indepen- 
dent, and mutually commute with one another, thereby 
constituting a complete set of integrals of motion. Thus, 
any hamiltonian that can be written in terms of these Ri 
likewise defines an exactly-solvable model. 

It is important to note that the set of integrals of mo- 
tion given in Eqs. © and (|10|) define a totally new set of 
exactly-solvable models, even though they were derived 
from the trigonometric family of RG models. That they 
are not simply the trigonometric family rewritten can be 
seen by focussing on Eq. (|10fl . which gives the form of the 
new Rj integrals of motion. They are in fact identical to 
those of the rational family of ESMs, except for the last 
two terms which are essential for ensuring the commuta- 
tion with the new bosonic integral of motion Rq. In this 
sense, the new class of models that we derived can be 
viewed as an extension of the rational family to include 
an extra boson degree of freedom. 

Note further that the operator that counts the total 

number of pairs, M = tfb + \ J2jm a jm a j*n> a ^ so com- 
mutes with all Ri and thus defines a conserved quantity. 



Before continuing our derivation of the properties of 
the exact solutions associated with this complete set of 
integrals of motion, we first write down some of the pos- 
sible hamiltonians that could be treated exactly in this 
way. One that is especially interesting is obtained di- 
rectly from the selected integral of motion Rq, 
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where V = wgWy and Ej — Vsj/2. In the pair repre- 
sentation of SU(2) or SU(1,1), this hamiltonian describes 
the interaction of fcrmionic or bosonic atom pairs with 
a diatomic molecule. In the two-level representation, it 
generalizes the Tavis-Cummings model to multi-atoms. 

Many other exactly-solvable models can also be con- 
structed in this way. For example, by taking a linear com- 
bination of the other integrals of motion Rj [j = 1, ..., L) 
with coefficients e,-, we obtain a hamiltonian of the form 
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Included are those that contain pairing interactions and 
those with level-dependent atom- molecule couplings. 

We now discuss how to rewrite the seniority-zero so- 
lutions for the trigonometric RG models to apply when 
one of its degrees of freedom has been replaced by an 
ideal boson in the infinite-Oo limit. Defining x a — 

y 2 cot ( e ° — 'To)' t ne Richardson-Gaudin equations (J5J 
that define the parameters e a , and thus the x a , become 
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The corresponding expressions for the eigenvalues © 
associated with the new integrals of motion I|9I10|1 are 
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while the seniority-zero eigenvectors take the form 
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We note here that each independent solution of the set on 
non-linear coupled equations l|13[l defines an eigenstate 
(|16fl that is common to the L + 1 integrals of motion 
(I9I1U|) and has eigenvalues I|14I15|I . 

The eigenvalues of the hamiltonian (|llf) . for example, 
can be obtained from the eigenvalues l|14|) of ro as 



E = VY,. 



(17) 



The eigenvalues of the hamiltonian l|12|) can similarly be 
obtained from the eigenvalues (|15fl of the R b operators. 

It is worth noting here that the solutions given in 1131 - 
| 16|l arc identical to those of the Tavis-Cummings model 
[3 for a single SU(2) spin and to those for the atom- 
molecule model of 9] for a single SU(1,1) bosonic level. 

Important observables in these models are the occupa- 
tion probabilities of the various degrees of freedom. They 
can be obtained from the integrals of motion using the 
Hellman-Feynman theorem, viz: 
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(12) from which we obtain for the occupations numbers 
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The derivatives of the x a are obtained by differentiating 
the RG equations (|T3)l with respect to G, which gives 
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We now turn to a specific application of these new 
exactly-solvable models as an illustration of the proce- 
dure to solve eq. I|13[) an d (|20() . We consider a mixture 
of bosonic atoms confined to a 3D isotropic trap cou- 
pled to a molecular two-particle bound state, a molecular 
dimer, and model it through the hamiltonian l|ll|l . This 
hamiltonian does not contain an atom-atom interaction, 
which could be included by using the more general hamil- 
tonian of l|12|) . In the hamiltonian we use, uj is the en- 
ergy of the molecular dimer above that of the Feshbach 
resonance and is the negative of the detuning parame- 
ter. Also, V is the atom-molecule interaction strength, 
e j = J (j ' = 0, 1, • ■ ■ ) are the single-atom energies in a 
3D isotropic trap, and Vtj = (j + l)(j + 2)/2 are the 
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as the interaction strengthens, but there is no phase tran- 
sition and the occupation of the lowest trap level is always 
negligible. A detailed study of the phase diagram that 
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FIG. 1: Occupation probabilities as a function of the inter- 
action V for a molecule energy uj = 10. The solid line corre- 
sponds to condensed atoms, the dashed line to non-condensed 
atoms, and the dotted line to condensed molecules. 



level degeneracies. The phase diagram and dynamics of 
this model have been studied in several recent works, e.g. 
[l8| . To make contact with ref. [jj| , our detuning param- 
eter is related to theirs by uj — S and our atom-molecule 
coupling is V = -^g • 

As noted before, the complete set of seniority-zero 
eigenstates arise from different solutions of the equations 
lfT3jl . For boson systems with negative couplings V, the 
parameters x a for the ground-state solution are real and 
positive. Excited states correspond to different partitions 
of the x a in intervals defined by the Sj . 

We have performed calculations for a system with 
M = 500 pairs and two values of the molecular energy, 
u> = 10 and uj — —10, as a function of the negative cou- 
pling V . The atom space was truncated to L — 50 har- 
monic oscillator shells. In Fig^ we show the occupation 
probabilities of the atomic condensate (solid line) and 
the atomic depletion (dashed line), and the fraction of 
molecules (dotted line) as a function of V for positive 
molecular energy (uj = 10). As can be seen, a quantum 
phase transition takes place at V ~ —0.45. Interestingly, 
the occupation of the atom condensed state is negligi- 
ble for V < —0.45 and the atomic fraction is distributed 
among all harmonic oscillator levels. A pure molecular 
state does not exist for any value of V. This quantum 
phase transition was recently studied using mean-field 
and renormalization group techniques and it was con- 
cluded that it lies in the Ising universality class [ijj . The 
exact solution offers a unique opportunity to unravel the 
critical properties around the transition point. 

In Figure |21 we show results for negative molecular en- 
ergy, uj = —10. The system is purely molecular for weak 
coupling. Molecules begin to decay to pair atomic states, 
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FIG. 2: Occupation probabilities as a function of the interac- 
tion V for a molecule energy ui = — 10. The solid line corre- 
sponds to condensed atoms, the dashed line to non-condensed 
atoms, and the dotted line to condensed molecules. 



emerges from the exact solutions will be given elsewhere. 

In closing, we have presented in this paper a new fam- 
ily of integrable models for atom-molecule systems. The 
models are exactly solvable for fermionic and bosonic 
atoms interacting with molecular dimers. There is a 
large freedom to select the parameters of the hamilto- 
nian, allowing for the description of quite general realis- 
tic systems. We have presented initial results for a mixed 
system of trapped bosonic atoms and molecular dimers. 
Application to systems of fermionic atoms and molecular 
dimers is of special interest due to the recent generation 
of an ultra cold molecular BEC from the conversion of 
40 Ka [b|] or 6 Li [lj| fermionic atoms. Such models can 
also be used to explore the BCS to BEC crossover, from 
a condensate dominated by Cooper pairs to a condensate 
dominated by molecular dimers. Finally, it is possible to 
use alternative realizations of these models to describe 
problems of importance in quantum optics and perhaps 
elsewhere. 

After this work has been submitted, we learn 
about a recent preprint on generalized integrable 
matter-radiation models [2(j. The extended the 
Jaynes-Cummings Hamiltonians treated there have non- 
hermitian atomic interactions. 

This work was supported by the Spanish DGI un- 
der grant #s BFM2003-05316-C02-01/02, by the US Na- 
tional Science Foundation under grant #s PHY-9970749 
and PHY-0140036, and by UBACYT X-204. Fruitful dis- 
cussions with L. L. Sanchez-Soto are acknowledged. 



[1] E. Jaynes and F. Cummings, Proc. IEEE 51, 89 (1963). [11] 

[2] G. Rempe, H. Walther, and N. Klein, Phys. Rev. Lett. [12] 

58, 353 (1987). [13] 

[3] M.G. Raizen et al, Phys. Rev. Lett. 63, 240 (1989). 
[4] M. Tavis, and F.W. Cummings, Phys. Rev. B 170, 379 [14] 

(1968) [15] 

[5] N.M. Bogoliubov, R.K. Bullough, and J. Timonen, J. 

Phys. A 29, 6305 (1996). [16] 

[6] A. Rybin, et al. J. Phys. A 31, 4705 (1998). 
[7] B. Buck, and C.V. Sukumar, Phys. Lett A 81, 132 (1981). [17] 

[8] LP. Vadeiko, G.P. Miroshnichenko, A.V. Rybin, and J. [18] 

Timonen, Phys. Rev. A 67, 053808 (2003). [19] 

[9] J. Links, H.-Q. Zhou, R.H. McKenzie, M.D. Gould, J. 

Phys. A 36, R63 (2003). [20] 

[10] J.N. Milstein, S.J.J.M.F. Kokkelmans, and M.J. Holland, 

Phys. Rev. A 66, 043604 (2002). 



K. M0lmer, Phys. Rev. Lett. 90, 110403 (2003). 
K. Xu, et. al., Phys. Rev. Lett. 91, 210402 (2003) 
M. Greiner, C. A. Regal, and D. S. Jin, Nature 426, 537 
(2003). 

S. Jochim, et. al, Science 302, 2101 (2003). 
J. Dukelsky, C. Esebbag, and P. Schuck, Phys. Rev. Lett. 
87, 066403 (2001). 

J. Dukelsky, S. Pittel and G. Sierra, submitted to Re- 
views of Modern Physics (2003). 
M. Gaudin, J. Phys. (Paris) 37 1087 (1976). 
M. Kostrun, and J. Javanainen, cond-mat/0308259 
L. Radzihovsky, and J. Park, cond-mat/0312237 M.W.J. 
Romans et al, cond-mat/0312446 
A. Kundu, quant-ph/0307102 



